Bound states in inhomogeneous magnetic field in graphene: a semiclassical approach 



00 

o 
o 

(N 

>>: 



A. Kormanyos,^ P. Rakyta,^ L. Oroszlany,^ and J. Cserti ^ 

^Department of Physics, Lancaster University, Lancaster, LAI 4YB, UK 
"^Department of Physics of Complex Systems, Eotvos University, 
H-1117 Budapest, Pdzmdny Peter setdny 1/A, Hungary 

We derive semiclassical quantization equations for graphene mono- and bilayer systems where 
the excitations are confined by the applied inhomogeneous magnetic field. The importance of a 
semiclassical phase, a consequence of the spinor nature of the excitations, is pointed out. The 
semiclassical eigenenergies show good agreement with the results of quantum mechanical calculations 
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I. INTRODUCTION 

Observation of massless Dirac fermion type excitations 
in recent experiments on graphene has generated huge 
interest both experimentally and theoreticallyii^i^. For 
reviews on graphene see Refs.^HHQIland a special issue 
in Ref. 9. 

The intensive theoretical and experimental work have 
lead to good understanding of the physical phenomena in 
the bulk of disorder-free graphene in homogeneous mag- 
netic fieldi^i^. Recently, the interest in inhomogeneous 
magnetic field setups has also appeared. Martino et alJ^ 
have demonstrated that massless Dirac electrons can be 
confined by inhomogeneous magnetic field and that a 
magnetic quantum dot can be formed in graphene, lev- 
els of which are tunable almost at will. The so-called 
"snake states" known from studiesi^i^ on two dimen- 
sional electron gas (2DEG) have also attracted inter- 
est and their properties have been discussed in graphene 
monolayer-i^ii^ and in carbon-nanotubesi^iiS. 

Semiclassical methods have helped our understanding 
of complicated physical phenomena enormously and be- 
come a standard tool of investigation. Not only they of- 
fer a simple, easy to grasp classical picture but in many 
cases they can also give quantitative predictions on ob- 
servables. Yet the first semiclassical study on graphene 
systems has only very recently appeared^. UUmo and 
Carmier derived and expression for the semiclassical 
Green's function in graphene and studied the "Berry- 
like" phase which appears in the semiclassical theory. 
The importance of "Berry-like" and " non-Berry- like" 
phases in the asymptotic theory of coupled partial differ- 
ential equations and their role in semiclassical quantiza- 
tion have previously been discussed in Refs. [2^l23ll2^l25l . 

In this work we study a graphene nanoribbon^^ 
in a non-uniform magnetic fiel d^^'^^i^^ as shown in 
Fig. [ija) and (b) and a circular magnetic dot in graphene 
monolayer^^ [see Fig. [TJc)]. We assume that the appHed 
perpendicular magnetic field of magnitude \Bz\ changes 
in step-function like manner at the interfaces of the mag- 
netic and non-magnetic regions and that it is strong 
enough so that the magnetic length Ib = \/h/e\Bz\ is 
much smaller than the characteristic spatial dimension 



of the graphene sample. We show that in this case the 
semiclassical quantization can predict and help to under- 
stand the main features of the the quantum spectra at 
the K pointii^i^i^ii^ of graphene mono- and bilayer. 




FIG. 1: In a) and b) the applied perpendicular magnetic field 
B is zero in the center region of width 2W. In the case of a), 
in the left and right regions the magnetic fields point in the 
same directions, while in the case of b), in the opposite direc- 
tions. The magnitude B of the magnetic field is the same in 
both left and right regions. In the case of Fig. [ijc) a circu- 
lar non-magnetic region of radius R is considered in graphene 
monolayer, whereas outside this region there is perpendicu- 
larly applied magnetic field of magnitude B. 



The article is organized in the following way. Firstly, in 
section we give a brief overview of the exact quantum 
mechanical treatment of graphene mono and bilayer. We 
also discuss some of the technical details of the quantum 
calculation regarding the system shown in FiglT] Next, 
in section IIIII we introduce our semiclassical formalism 
and wherever possible we give a unified description for 
graphene mono and bilayer. In section IIVI we present 
the results of the semiclassical quantization for graphene 
nanoribbons in inhomogeneous magnetic field and com- 
pare it to tight binding (TB) and exact quantum calcu- 
lations. In the next, section|V]we apply the semiclassical 
formalism to a magnetic dot in graphene monolayer. Fi- 
nally, in section IVll we arrive to our conclusions. 
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II. MONOLAYER AND BILAYER GRAPHENE: 
QUANTUM MECHANICAL TREATMENT 

In the simplest approximation the Dirac Hamiltonian 
describing the low-energy excitations at the K point of 
the Brillouin zone in monolayer (bilayer) graphene reads 



(1) 



Here /? 1 (2) for monolayer (bilayer) and gi — vp — 
V3/2 ail/ft is given by the hopping parameter ti between 
the nearest neighbours in monolayer graphene (a = 0.246 
nm is the lattice constant in the honeycomb lattice). 
Moreover, 52 = — 1/2to* and the mass term is given by 
TO* = t2/2(5i)^, t2 being the interlayer hopping between 
A — B sites of bilayer—. The operators {Ttx,ny) are de- 
fined by TT = [tTx, T^y) = p + eA, where p = |^ is the 
canonical momentum operator and the vector potential 
A is related to the magnetic field through B = rotA. 
Due to the chiral symmetry azHpaz = —Hfj, where <Tz 
is a Pauli matrix, it is enough to consider the positive 
eigenvalues of the Hamiltonian Hfj . 

In section IIVI we study a graphene nanoribbon'^^ of 
width L ^ Ib [see Figs. [2a),(b)]. In the central re- 
gion |x| < W the magnetic field is zero, while for 
|a;| > a non-zero perpendicular magnetic field is ap- 
plied. We assume a step function-like change of the 
magnetic field at x — W and use the vector potential 
A(r) = {0,Ay{x),0)'^ to preserve the translation invari- 
ance in the y direction. Other details of the quantum 
calculation can be found in Refs. [T^flslflGl . 

The magnetic dot system is shown in Fig. [ijc). It 
consists of a graphene monolayer in homogeneous per- 
pendicular magnetic field with a circular inclosure where 
the magnetic field is zero. The circular symmetry of the 
setup suggests that one should choose the vector poten- 
tial in the symmetric gauge: 



A(r) 




(2) 



where r = (r cos (p, r sin 1^9) is in polar coordinates. One 
can show that with this choice the Schrodinger equation 
Hjs"^ p — E'^ 13 becomes separable in r and In the case 
of graphene monolayer, requiring that the wave function 
be normalizable and continuous a.t r — R leads to a secu- 
lar equation, solutions of which are the quantum eigenen- 
ergies (see Eq. (21) in Ref. Ila). 



III. SEMICLASSICAL FORMALISM FOR 
GRAPHENE 

We now give a brief account of our semiclassical for- 
malism. Our discussion goes along the lines of the 



Refs. [mili from where we have also borrowed some 
of the notations. 

We seek the solutions of the Schrodinger equation 
Hp^ p = p in the following forro^^: 



fc>0 



(3) 



where a^(r) are spinors and Spiv) is the classical ac- 
tion. Performing the unitary transformation p — > 
e-^s^("-)^-^, Hp e-^^'^("")i/;3e*^f»("") the Schrodinger 
equation can be rewritten as 

-E 



gp{Tlx~myf \ / p 
gp{Il.x + tny)P -E jr^^''^ 

-a?(r) + .. 



(4) 



= 0. 



Here H, - P.+HI, where H^ = px + eA,{r), = 
and similarly for Hj,. The WKB strategy^ is to satisfy 
Eq. ^ separately order by order in h. 
At 0{hP) order we obtain 



-E gpiUl fflO)'^ 



gp{U'i + ^m0 



'E 



a(?(r)=0. (5) 



This classical Hamiltonian can be diagonalizcd with 
eigenvalues 



7^±(p,r) = ±5M(n",(r))2 + (H°(r))Y 



(6) 



and eigenvectors V^(p,r). What we have found is that 
the 0{h?) order equation is in fact equivalent to a pair 
of classical Hamilton- Jacobi equations: 



E-n^p 



dv 



,r =0. 



(7) 



The solution of Eq. (l7])-when it exists- can be found e.g. 
by the method of characteristics^i. 

For £' 7^ the eigenvectors of the classical Hamiltonian 
given in Eq. ^ are 



V2 V 



(8) 



(here cj) is the phase of H^ — zH^) but the eigenspinor ap'^ 
can be more generally written as ap'* = A^{v)e^'^f 
where A^{r) is a real amplitude and 7|^(r) is a phase. 
Equations for ^^(r) and 7^(r) can be obtained from 
the 0{h^) order of Eq. (g]). One can show that the 0{h^) 
order equation can be written in the following form^iSl!^ 



{4'^)^Mp4'^ = 0. 



(9) 



Using the notation cr = {ax^cTy), <yx,y,z being the Pauli 
matrices, the operator Mi for graphene monolayer is 
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Ml = crp, while for bilayer it reads M2 = rh + rh^ , where 

7n = ^Tp(no + za,nO). 

The imaginary part of Eq. ([9]) expresses the con- 
servation of probability since it can be cast into the 
form of a continuity equation divj^ — 0. Here — 

Im('I'^'^|v^|^'^'*) is the probability current carried by 
the semiclassical wavefunction 5*^'^ — a^'^e'^^f ( = 

jr[Hf3,r] is the velocity operator). Similarly to the case 
of quantum systems described by scalar Schrodinger 
equation'^ «, this continuity equation determines (r) . 

The real part of Eq. ^ allows to calculate the phase 
7^ (r) . The equation determining 7^ (r) reads 



rf7^(r) 




gng(r) 

dy 



2 



where 



(10) 

±l/m*, and we denoted by 



Bz{v) the component of the applied magnetic field which 
is perpendicular to the graphene sheet. The second 

equahty in Eq. ^ follows from ^-^^ = ■ 

The Hamiltonian given in Eq. ([1]) yields a gapless spec- 
trum for bilayer graphene. Theoretical and experimental 
studies of bilayer graphen o^^'^^ have shown that an elec- 
tron density dependent gap can exists between the oth- 
erwise degenerate valence and conductance bends (de- 
scribed by Ti.2 and Ti.2 respectively, in our semiclassi- 
cal formalism). Assuming that the gap is spatially con- 
stant, one can take it into account by considering the 
Hamiltonianii H2,a — H2 + {/S./2)az- The new term 
(A/2)(7z affects only the 0{1tP) calculations, while the 
operator M2 in Eq. ^ remains the same. Consequently, 
Eq ([6]) is modified to 



H2±(p,r) =± 



1 



2m* 



A2 



(11) 



where A = A/to* and the right hand side of Eq. pn]) for 
/3 = 2 is multiplied by 77 = 1 — 



1 

It has been shown in Refs. |2^i24i25| that for iV- 
dimensional integrable systems where the particles have 
an internal, e.g. spin or electron-hole degree of freedom, 
one can derive a generalization of the EBK'^'^ quantiza- 
tion of scalar systems. In general, the quantization con- 
ditions read 



p dr 



In (n^ 
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(12) 



Here Tj,i — \...N are the irreducible loops on the A^- 
torus in the phase space, rijS are positive integers, /ZjS 
are the Maslov indices'^- counting the number of caustic 
points along Tj and finally a^s measure the change of the 
phase of the spinor part of the wavefunction as the sys- 
tem goes around a loop Tj . The systems we are consider- 
ing (see Fig. [T]) are more simple in that the Schrodinger 
equation is separable if the vector potential A(r) is cho- 
sen in an appropriate gauge, which takes into account 



the symmetry of the setup [i.e. translational symmetry 
in the case of Fig. [Ija), (b) and rotational in the case 
of Fig. [He)]. The magnetic field in Eq. ^ will de- 
pend on only one of the (generalized) coordinates, let us 
denote this coordinate by xi, the conjugate momentum 
by pi, and the other coordinate (conjugate momentum) 
by X2 {p2)- It turns out that due to the symmetry of 
the system 7(r) will also depend only on xi^^ . Therefore 
one of the two quantization conditions, involving the co- 
ordinate X2 and the conjugate momentum p2, is exactly 
the same as it would be for a scalar wavefunction [this 
corresponds to a2 = in Eq. (jl2p ]. In the quantization 
condition involving pi and xi however, the phase ai is 
in general not zero but is determined by Eq. (fTO|) : 



ai 



Bz{xi{t))dt. 



(13) 



For systems with piece- wise constant magnetic field pro- 
files such as those shown in Fig. [U the calculation of 7^ 

simplifies to 7^ = f X); Here T/ is the time 

the particle spends during one full period of its classi- 
cal motion in the Zth region where the strength of the 
perpendicular component of the magnetic field is given 
by Bz,i- In the semiclassical picture 7(r) changes only 
when the particle, during the course of its classical mo- 
tion, passes through non-zero magnetic field regions and 
this phase change of the wave function needs to be taken 
into account in the semiclassical quantization. 



IV. BOUND STATES IN GRAPHENE 
NANORIBBONS 



We now apply the presented semiclassical formalism to 
determine the energy of the bound states in inhomoge- 
ncous magnetic field setups in graphene nanoribbons, see 
Figs[T](a) and (b). Throughout the rest of the paper we 
will only consider corresponding to positive energies. 

would describe negative energies, these however do 
not need to be considered separately due to the chiral 
symmetry of the Hamiltonian, as explained in section HTl 
Using the Landau gauge A = (0, Ay{x), 0)-^ the trans- 
lation invariance of the system in the y direction is pre- 
served and therefore the solution of the Hamilton- Jacobi 
equation Eq. ([7]) can be sought as ^^(r) = Sp{x) +Pyy, 
where py = const. Since the classical motion in the y 
direction is not bounded, Py = hky is not quantized, it 
appears as a continuous parameter in our calculations. 
In contrast, the motion in the x direction is bounded due 
to the X dependent magnetic field Bz{x). Therefore the 
quantization condition reads 



^ f Pi3{x)dx + J/3 



27r(n-f 1/2). 



(14) 



(Note that the Maslov index is /i = 2.) It is useful to in- 
troduce at this point the following dimensionless param- 
eters: the width of the non-magnetic region w = W/Ib, 
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the guiding center coordinate X = kyls, both in units of 
Ib (which is defined in Sec. |TT]). Throughout this paper 
we will use w = 2.2. 

We start our discussion with the magnetic waveguide 
configuration shown in Fig. [TJa) in graphene monolayer. 
Introducing the dimensionless energy Emi = Els/hvp, 
one finds that for \X\ < Emi there is one turning point 
in each of the left and right magnetic regions. A simple 
calculation gives 71 = tt and writing out explicitly the 
result of the action integral from Eq. it follows that 

4.Kmiw + 7rEl^i^2nn, ^==1,2,.... (15) 

Here we h ave introdu ced the dimensionless wave number 
Kmi — \J E^i — X"^ and note that the phase change of 
the wave function due to 7(0;) cancelled the phase con- 
tribution coming from the Maslov index. Furthermore, 
if X > Emi {X < —Emi) there are two turning points 
in the left (right) magnetic regions. One finds that also 
for this case 71 = tt which again cancels the contribution 
from the Maslov index, thus for \X\ > Emi the semiclas- 
sical quantization yields 

£'„ = V2^, n = l,2,... (16) 

i.e. the energies are independent of X (and hence 
ky). This is the same as the exact quantum and the 
semiclassical2i results for the relativistic Landau levels 
(LLs) in homogeneous magnetic field. [From the exact 
quantum calculations^^' it is known that a zero energy 
state also exists in this system. Formally, from Eq. pB]) 
one can obtain a zero-energy state by assuming that 
n = is admissible. However, Eq. ([5]) and hence Eq. (fTUl) 
are only valid for i? 7^ 0. Therefore we exclude n = 0.] 

Comparison of the semiclassical eigenvalues with the 
results of exact quantum calculations are shown in Fi g. El 
(For details of the quantum calculation see e.g. Ref. KLg .) 
The agreement between the quantum and semiclassical 
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FIG. 2: Results of exact quantum calculations (solid lines) 
and the semiclassical approximation given by Eqs. (|15|) and 
P6[) (circles) as a function of ky (in units oi Ib) for graphene 
monolayer. The dashed lines indicate \X\ — E,ni (see text). 

calculations is in general very good, especially for higher 



energies. For lower energies and |^| ^ E^i one can ob- 
serve quantum states which are almost dispersionless and 
their energy is very close to the non-zero energy LLs in 
graphene monolayer. Semiclassically, these states are de- 
scribed by Eq. . Although the zero energy state of the 
spectrum^^ can not be accounted for by our semiclassics 
an expression for the gap between the zero and the first 
non-zero energy states can be easily obtained by putting 
n = 1 and AT = in Eq. (jlSp and it gives a rather accurate 
prediction as it can be seen in Fig.O The presented semi- 
classical method can not describe those quantum states 
which correspond to \X\ « Emi [see the dashed line in 
Fig. [2] i.e. when one of the the turning points is in the 
area of rapid spatial variation of the magnetic field. 

For comparison, we have also calculated the quanti- 
zation condition for graphene bilayer using the classical 
Hamiltonian given in Eq. (jll|) and the general quantiza- 
tion condition shown in Eq. (HH). For Eu > X'^/2, where 
Ehi — uic — -^^^ being the cyclotron frequency, it 
reads: 

2Km w + nEu =^(j^-l^ ^ 71 = 2, 3, ... . (17) 

Here Ku = \/ 2Eu — X'^ and we have taken into account 
that in this case 72 = 27r. For rj = 1 [where rj has been de- 
fined after Eq. (|lip ] this result is very similar to what one 
would obtain for a 2DEG, the only difference being that 
for 2DEG one would have -1-1/2 on the right hand side 
of Eq. fTT]) . This similarity is a consequence of having 
a parabolic dispersion relation £^(k) close to the Fermi 
energy in both a 2DEG and graphene bilayer systems. 
We let the integer quantum number n to run from n = 2 
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FIG. 3: Results of exact quantum calculations (solid lines) 
and the semiclassical approximation given by Eqs. (|17|l and 
(|18p (circles) as a function of ky (in units of Zs) for graphene 
bilayer. The dashed lines indicate — Eu (see text). 

in Eq. (fT7|) for the following considerations: from Ref. [21I 
we know that in a more simple case of homogeneous mag- 
netic field the four-fold degenerateS^ LL of the quantum 
calculationgii at Eu = (corresponding to n = 0, 1) can 
not be correctly described semiclassicaly but for LLs hav- 
ing Eu > the agreement between the semiclassical and 
quantum results is qualitatively very good. Similarly, we 
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expect that in our case the semiclassical approximation 
should only work for n > 2. We have found that this 
is indeed the case, see Fig. |3] where the solid lines show 
bands obtained by TB calculations and the circles are cal- 
culated using Eq. (fTT)) for n > 2 [we have taken rj — 1]. 
The Ebi > energy bands for Em > X'^/2 are remarkably 
well described by Eq. ([TT]) [note however that like in the 
homogeneous magnetic field case, there is a four-fold de- 
generate state at Em = ]. For Em < X^/2 the bands of 
TB calculations again become almost dispersionless and 
level off very close to the LLs of bilayer graphene in ho- 
mogeneous field--l. The semiclassical expression for the 
energy levels in this regime of X is 



Eh, 



(n-1/2), n = 2,3. 



(18) 



which is again a good approximation of the quantum re- 
sult. Our semiclassics can not correctly account for states 
having X'^ /2 sa Em i.e. when one of the turning points is 
in rapidly varying magnetic field region. 

We now turn to the semiclassical study of the system 
depicted in Fig.[Tfb), where the magnetic field is reversed 
in one of the regions. It has been shown in Refs. 
that peculiar type of current carrying quantum states 
called snake states exist close to the K point of graphene 
for this magnetic field configuration. These states can 
also be described by the Dirac Hamiltonian and are there- 
fore amenable to semiclassical treatment. 

We start the discussion with the graphene monolayer 
case. There are no turning points and hence no states if 
—X > Erni- For \X\ < Emi there is one turning point in 
each of the non-zero magnetic field regions. In contrast 
to the symmetric magnetic field configuration [Fig. [IJa)] 
we find that for one full period of motion 71 = 0: the 
contributions of the two magnetic field regions, point- 
ing in the opposite direction, cancel. Using Eq. the 
quantization condition is: 



7r(n + 1/2) 





" X ' 




^arcsin 




•5) 




Eml 





n = 0,1,2,... 

(19) 

Note, that unlike in the case of Eq. (|15p here a solution 
for n = also exists. Moreover, for X > E^i there are 
two turning points both in the left (we denote them by 

X2 leads to the same 
/2n, while using xf. 
The difi'er- 



1) and in the right (denoted by ) magnetic 

region. The quantization using 
result as in Eq. i.e. En,L = 
X2 gives the sequence En,R = \j2{n + 1 
ence between En^L and En.R is due to the fact that the 
sign of the phase contribution from 7(0;) depends on the 
direction of the magnetic field, i.e. it is -|-7r when xf, 
x\ is used in the calculations and — tt when ccf , x^ is 
used [see Eq. (ITUl) ]. From these considerations it follows 
that if 71 = 0, 1, 2, . . . as we assumed in Eq. ([T^. the two 
sequence En^L, En^R give two- fold degenerate dispersion- 



less states at E^i ~ \f2, VZ, ^/6 . . . and a nondegener- 
ate one at Emi = 0. We have to exclude, however, the 
Eml ~ solution, see the discussion below Eq. fTB]) . 

The results of quantum and of semiclassical calcula- 
tions are shown in Fig. [D (For details of the quantum 




FIG. 4: Results of exact quantum calculations (solid lines) 
and the semiclassical approximation given by Eq. (I19|l (cir- 
cles) as a function ky (in units of Ib) for graphene monolayer. 
The dashed line indicate X = Emi- (see text). 

calculation see e.g. Ref.y.) As one can see the agree- 
ment is again very good for Emi ^ 0.4 except when 
X « Eml, see the discussion in the previous example, 
i.e. when the magnetic fields point in the same direction. 
The two-fold degenerate dispersionless quantum states 
can also be observed for X > Emi ■ States having energies 
< Eml ^ 0.25 can not be described by our semiclassics 
and note that in the lowest energy band corresponding to 
n = in Eq. (fT9|) nonphysical solutions also appear along 
with the genuine ones for 0.25 < Emi ^ 0.4. This clearly 
indicates the limits of applicability of our approach, i.e. 
it docs not work for energies close to the Dirac point. 

We end our discussion of the bound states in mono 
and bilayer graphene nanoribbons with the bilayer sys- 
tem corresponding to the previous, monolayer example, 
e.g. for the magnetic field setup of Fig. [Ijb). For 

-X > ^/2E M there are no turning points and hence no 
states. The quantization condition for X^/2 < Em can 
simply be obtained from Eq. (jl9p by changing Emi ^ 
\/2Em and Kmi Km \Km is defined after Eq. (fTT])]. 
Finally, for X"^ /2 > Em our semiclassics predicts a se- 
quence of doubly-degenerate, dispersionless energy lev- 
els at Em — [n — 1/2), n = 2, 3, . . . , in a similar fash- 
ion as in the monolayer case. As one can see in Fig. [S] 
for Em ^ 0.8 the semiclassical approximation captures 
all the main features of the TB calculations quite well, 
apart from the region where X^/2 w Em for X > [see 
the discussion below Eq. (HH])]. Dispersive states corre- 
sponding to n = 0, 1 can also be described semiclassi- 
cally if Em > 0.8 [see the lowest two bands in FiglS], but 
for smaller energies nonphysical solutions along with the 
genuine ones do appear and for Em ^ 0.1 i.e. very close 
to the Dirac point no quantum states can be described 
with the presented semiclassical approach. 
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FIG. 5: Results of TB calculations (solid lines) and a semi- 
classical approximation (circles) for graphene bilayer as a 
function of ky (in units of Ib)- The semiclassical approxima- 
tion can be obtained from Eq. p9p by a transformation de- 
scribed in the main text. The dashed line indicate /2 — Eu 
(see text). 



V. BOUND STATES OF A MAGNETIC 
QUANTUM DOT 

Our last example is the magnetic dot in graphene 
monolayer discussed in Ref. [H and shown in Fig. [UJc). 
We assume that the magnetic field is zero in a circular 
region of radius R while outside this region a constant 
perpendicular field is applied. 

Working in polar coordinates r and ip^ since the vector 
potential of Eq. ([2]) preserves the circular symmetry of 
the problem, the Lp coordinate is cyclic and therefore one 
can seek the solution of Eq. ([7]) as ^(r) — Sr{r) + PipP 
where p^p — const. Moreover, in ip coordinate the motion 
is free rotation and hence the quantization condition for 
p^ is simply 



Pipdkp — 2h'Km, 



0,±1,±2. 



(20) 



whence it is clear that the pip quantization reads p^p = 
hm. 

The second quantization condition can generally be 
written as j- f pr{r,m)dr -I- 71 = 27r(n + 1/2) because 

the Maslov index is = 2 and pr = We now 

2 

introduce the dimensionless variable ^ = ttpt and the pa- 

rameters S — = i?^/2, m — m — S. One can see that S 
is basically the missing magnetic flux that can be associ- 
ated with the dot. The phase accumulated between two 
points ^1, ^2 inside the dot (where Bz{r) = 0) is given by 



S?=\iu^2)^lj ^ (21) 



while between two points in the non-zero magnetic field 



region by 



-^2 + 2(i^2 „ . 



(22) 

[the h factor in the above expressions appears because in 
the calculation we already took into account the quan- 
tization of see Eq. (pil)) ]. The integrals in Eqs. (PT|) 
and ((22)) can be analytically calculated but the resulting 
expressions are too lengthy to be recorded here. 

As a next step to obtain a semiclassical quantization 
rule we proceed with the analysis of the classical dynam- 
ics along the lines of Ref. |2^ Calculating the radial ve- 
locity Vr{r) = ^^^Pr in the magnetic region one finds 
that for i?2 > 2m there can be two turning points in the 
radial motion which we denote by , ^(^. In terms of the 



dimensionless parameters Rc = E and X = y R"^ — 2rh 
(the radius of the classical cyclotron motion and the guid- 
ing center coordinate respectively, both in the units oUb) 
the turning points can be written as = ^{Rc ± X)"^. 
With regard to S, there are then two possible cases. 

i) The first case is when S < Co" > (or equivalently, 

2 

S < and therefore the radial motion is confined en- 
tirely to the magnetic region. Calculation of 71 gives 
a phase +tt which cancels the phase contribution from 
the Maslov index. Hence the quantization condition is 
25,^^"(Co",Ccj")/ft = 2n7r where S^^^i^^^^o) is calcu- 
lated using Eq. (p2)) . Explicitly, the energy levels En^m 
are given by 



E„ 



m. 



(23) 



This result is very similar to what one would obtain from 
exact quantum calculations for a homogeneous magnetic 
field, where the relativistic Landau levels are given by 
El"^ = ^J2n+ \m\ + m. Note however that in Eq. ^ 
instead of the integer quantum number m the non-integer 
rh ^ m — S appears. From Eq. (|23[) it is clear that for 
m < the energy bands are at En = ■\/2n, n = 1,2,... 
and they do not depend on 5 and m. (As in the nanorib- 
bon case, we exclude n — Q because that would give 
En — 0.) These 5 independent sections are readily ob- 
servable in the energy bands corresponding to m = — 2 
in Fig. El^a). On the other hand, for to > 0, (5 ^ m an 
approximately linear dependence on 5 of the bands corre- 
sponding to different to is predicted by Eq. (j23p and this 
can also be observed, see Fig. [TJa). From classical point 
of view in the parameter range to < (which implies 
X — Rc > R) the classical orbits are such that they do 
not encircle the zero magnetic field region [see Fig.|6l[b)]. 
They are just like orbits in homogeneous magnetic field 
and this helps to understand why the quantum states cor- 
responding to the same parameter range are reminiscent 
of dispcrsionless Landau levels. Conversely, for to > 
(which implies Rc — X > R) the classical orbits do en- 
circle the zero magnetic field region and therefore the 
energy of the corresponding quantum states depend on 
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FIG. 6: a) results of exact quantum calculations (solid lines) 
for the m = —2 energy bands as a function of the missing flux 
S. The results of the semiclassical quantization (*) obtained 
from Eq. (dSll for £ < i^, m < 0. For E > circles 
(o) show the semiclassical results calculated using Eqs. (|24[) 
and (pSj) . The dashed line shows the E — function which 
separates the cases i) and ii) detailed in the main text, (b) 
shows a cartoon of a classical orbit in the parameter range 
E<i^,m<0. 

the missing flux 5 [see the cartoon shown in Fig. [Tljb) for 
illustration of the classical orbits] . 

ii) The second case is when < S < which hap- 
pens if (5 > The classical motion is no longer con- 
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FIG. 7: (a) results of exact quantum calculations (solid lines) 
for the m = 2 energy bands as a function of the missing flux 
S. The results of the semiclassical quantization (*) obtained 
from Eq. ((23ll for i5 < m > 0. For i5 > ^ circles 

(o) show the semiclassical results calculated using Eqs. (|24|) 
and (1251) . The dashed line shows the E — -fe= function which 
separates the cases i) and ii) detailed in the main text, (b) 
shows a cartoon of a classical orbit in the parameter range 
E<^,fh>0. 

fined to the magnetic region but also enters the nonmag- 
netic dot. The turning point in the nonmagnetic region 

2 

is at ^0 — < S. Since the Maslov index is fj, — 2, the 
quantization condition can be written as 

I [5f="(eo,<5) + 5,^^°(<5,eo+)] +71 = 27r(n+l/2). (24) 



Here S^=°{£.o,S) and S^^°{S,£_+) can be calculated us- 
ing Eq. ([2T|) and Eq. ([22]) respectively, but the resulting 
expressions are again too lengthy to be presented here. 
Moreover, from Eq. we find that 



71 = - + arcsin . = ■ (25) 

^ \^Je^{& -2m) J 

One can see that here in general 71 ^ tt and therefore it 
does not cancel the contribution of the Maslov index. 

As one can observe the overall agreement of the exact 
quantum and of the semiclassical calculations shown in 
Figs El [3 is good, especially for higher energies. Accord- 
ing to the exact quantum calculations^^, there is also a 
zero energy state but this can not be described by our 
semiclassics. 



VI. CONCLUSIONS 

To conclude, using semiclassical quantization we have 
studied the spectrum of bound states in inhomogeneous 
magnetic field setups in graphene mono and bilayer. We 
have found that a semiclassical quantization which takes 
into account a "Berry-like" phase can indeed explain all 
the main features of the exact quantum or numerical TB 
calculations. In particular, we have studied graphene 
mono-and bilayer nanoribbons in magnetic waveguide 
configuration and also in a configuration when snake 
states can exist. Besides, we discussed the magnetic dot 
system in graphene monolayer. For the considered step- 
wise constant magnetic field profile we have derived semi- 
classical quantization equations. In the case of graphene 
monolayer, we have compared the resulting semiclassical 
eigenenergies to quantum mechanical ones obtained from 
the corresponding Dirac equation. For graphene bilayer 
the results of the semiclassical quantization and numeri- 
cal TB calculations have been compared. In all the cases 
a good agreement has been found except for energies very 
close to the Dirac point. We have shown that the main 
features of the spectrum depend on whether the classical 
guiding center coordinate is in the non-magnetic or in 
the magnetic field region. 

Assuming homogeneous magnetic field, the energy of 
the Landau levels in semiclassical approximation has 
been calculated in Ref. Our work can be consid- 
ered as a generalization of these calculations to a class 
of non-homogeneous magnetic field setups, where due to 
the symmetry of the system, the Berry-like phase ap- 
pearing in the semiclassical theory affects only one of the 
quantization conditions. 
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